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The rotationally-covariant renormalization group equations of motion for the density wave ampli- 
tudes in the phase field crystal model are shown to follow from a dynamical equation driven by an 
effective free energy density that we derive. We show that this free energy can be written purely as a 
function of the strain tensor and thence derive the corresponding equations governing the nonlinear 
elastic response. 
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Multiple scale approaches to materials pattern forma- 
tion are essential to account for the variety of material 
structures that emerge on different scales, and the de- 
pendence of materials properties on structures at more 
than one scale However, a coarse-grained description 
describing the pattern dynamics is not in itself sufficient: 
for example, in fracture mechanics, one would like to be 
able to predict the response of real, heterogeneous ma- 
terials to deformations, and to couple the resulting elas- 
ticity theory to discrete, atomistic descriptions at small 
scales. In short, the challenge of materials modeling re- 
quires not just a coarse-grained description of the un- 
derlying atomistic dynamics, but a multi-resolution ap- 
proach that can locally adapt to capture the appropriate 
fast space and time scales. 

Recently, we presented a proof-of-concept realization 
of such a calculation 0, 0, 0], using a minimal model 
of a crystalline material called the phase field crystal 
(PFC) modelfl, i, Q.The PFC model evolves the local 
density field through a conservation law governed by a 
free energy functional that penalizes departures from a 
perfectly periodic ground state [1]. The resulting density 
field retains the crystallographic and elastic properties 
of the material, in contrast to conventional phase field 
model descriptions of materials, and appears to give a 
realistic account of a variety of phenomena, including 
multicrystalline solidification Q, elasticity@,[y,Q, defect 
dynamics [H, [l(| EH, epitaxial growth, as well as crack 
and fracture dynamics^. Our proof-of-concept calcula- 
tion consisted of two distinct steps: (i) a renormalization- 
group calculation for the rotationally-covariant equation 
of motion for the amplitudes of density waves correspond- 
ing to the lattice periodicity, and (ii) the numerical solu- 
tion of the resulting equations using adaptive mesh refine- 
ment techniques 4]. In two dimensions, this calculation 
was three orders of magnitude faster than straightforward 
integration of the PFC equations. 

In this Rapid Communication we show that an addi- 
tional level of coarse-graining can be performed on the 
complex amplitude representation, leading to a theory 
of nonlinear elasticity valid at system-wide or macro- 
scopic scales. Our construction is of interest method- 
ologically: instead of postulating a form for the strain 



energy based on symmetry and phenomenology [1^, Il3j |. 
we derive it from a more microscopic model — in this case, 
the PFC model. Our work provides a theoretical connec- 
tion to microscopic structures and opens up the possi- 
bility of deriving other forms of strain energy by similar 
approaches. Finally, in varying our level of description 
from the PFC model to the complex amplitude represen- 
tation to the strain energy formulation, we obtain the 
basis for a multiple scale approach to materials proper- 
ties, accessing the realm of continuum mechanics from 
a density functional theory at atomic scales. This ap- 
proach could be applicable to interpreting recent exper- 
iments on atomic force microscope nanoindentation of 
graphenefbfl, IT1|. complimenting theoretical approaches 
using tight-binding methods 16]. 

The Model:- The phase field crystal (PFC) model is de- 
fined by the free energy density d, Q, 



(1) 



where p(x, t) is the phase field, or the order parameter. 
The dynamics is conservative and dissipative, given by 



dp 
dt 
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(2) 



where F = j d 2 xf(x) is the total free energy of the sys- 
tem, taken here to be two-dimensional. Thermal noise 
is ignored in our discussion here, but could be included 
if desired. Generally, it will not be important at sys- 
tem scales. There are three phases in this model, namely 
uniform, stripe and triangular. The ground state in the 
triangular phase can be written in the single-mode ap- 
proximation by 



p tn {x) = ( 



Po, 



(3) 



where po is the average density, A is the constant ampli- 
tude and 



— X ~2 y (4) 
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are the three lattice vectors. Goldenfeld et al. showed [2, 
3] that instead of using p(x, t) as the dynamical variable, 
it is more efficient to generalize Eq. ([3]) and promote 
the constant amplitudes A to 3 slowly varying complex 
amplitudes Aj (x, t) and treat the Aj (x, t) as dynamical 
variables. Using renormalization group techniques, they 
showed that the rotationally-covariant equation of mo- 
tion for Ai(x, t) is given by[2j,|3| 



dAi 



(l-ii)(r-i?)Ai-6po^ 
-3A 1 (\A 1 \ 2 + 2\A 2 \ 2 + 2\A 3 \ 2 ) + ---, (5) 



where r = — r — 3/?q and L\ = V 2 + 2ik 1 • V is a rota- 
tionally covariant operator, and small nonlinear gradient 
terms [3| are not written out explicitly. The correspond- 
ing equations for A 2 ^(x, t) can be written down from the 
appropriate permutations. 

Derivation of the amplitude equations without using 
renormalization group:- The first attempt to derive the 
free energy in this representation is to note that Eq. |(SJ), 
ignoring the nonlinear gradient terms, can be written as 



dt 



8F dyn [Aj{x,t)} 
SA* 



(6) 



where the free energy Fd yn is given by the free energy 
density 

hyn = -^^(l-^r-i^. + a^ \A 3 \ 2 \A k \ 2 



3 = 1 



3,k=l 



~ 2 E l^l 4 + GMAiA 2 A 3 + A{A*A* 3 ), (7) 
3=1 

The dynamics given by Eq. ^ is purely dissipative, as 
opposed to the density-conserving dynamics in the orig- 
inal PFC equation, Eq. ©. 

To resolve this conundrum, we note that mass conser- 
vation should only be exhibited in the dynamical equa- 
tion of motion and not be represented in the equilibrium 
free energy. Thus, the correct way to derive the free en- 
ergy in the complex amplitude representation is to derive 
it from the original PFC free energy, Eq. |T| . The easiest 
way to do that is to substitute the ansatz, Eq. into 
the free energy, Eq. ([I]). The first term of Eq. (JTJ can 
then be computed by using the identity 

(1 + V 2 )p = J2 (' "•' %•>, + cc) + po, (8) 

3=1 

where c.c. stands for complex conjugate. By performing 
an integration by parts, we find 



f(l + V 2 ) 2 p = i[(l + V>] a 

= lt A 3 L 3 A 3> 
3=1 



(9) 
(10) 



where in the last line constants and terms with the 
rapidly oscillating factor exp (ikj ■ x) are neglected. We 
can neglect the oscillating terms because the complex 
amplitudes, A,, are slowly varying on that scale, so 
the terms cancel themselves upon integration over space. 
Other terms in the free energy can be transformed in a 
similar fashion. The resulting free energy is given by 



f amp = -J2A*(T-L*)A j + 3j2\ A 3\ 2 \M 2 

3=1 3,2=1 
3 3 

"2 E \ A i I 4 + ®Po(AiA 2 A 3 + A\A* 2 At){ll) 
3=1 

Note that this free energy is different from Eq. Q. The 
1 — Lj operator in the first term in Eq. ([7]) is absent 
here. This is to be expected because this operator arises 
from the conservative Laplacian in the dynamical equa- 
tion, and according to our discussion above, it should not 
appear in the free energy. The transformation of the dy- 
namical equation, Eq. p]). into the complex amplitude 
representation can be performed by observing that for 
any function, f(x), the identity, 



dPxftLj - l)/(x)]. 



0. 



holds, a counterpart to the identity, 

d 2 xV 2 g(x) = 0, 



(12) 



(13) 



for any function g{x). In fact, if we define g{x) = 

f(x)e^- s , Eq. (USD implies Eq. (JTSJ). This shows that 
when we make the change of variables from the density, p, 
to the complex amplitudes, Aj, the Laplacian in the con- 
servative dynamical equation has also to be transformed 
to Lj — 1. We thus arrive at the equation of motion 



dAj 
dt 



SF 

( j 1 \ ur amp 



(14) 



which, when written out explicitly, is, 



^ = (1 
dt 



L x )[{T-L\)Ai-&p Q AlAt, 
-3^ 1 (|^ 1 | 2 + 2|A 2 | 2 + 2|^ 3 | 2 )] (15) 



with appropriate permutations for A 2t3 (x,t). By con- 
struction, these equations conserve the density of the sys- 
tem. 

Eq. (fT5|) is exactly the same as Eq. {!J| with all the 
nonlinear gradient terms included. This derivation shows 
that the inclusion of the nonlinear gradient terms is cru- 
cial for density conservation, and that all those terms can 
actually be written in the condensed form of Eq. (|15p . 
Note, however, that our derivation does not explicitly 
use a renormalization group argument, but follows from 
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the integral identity, Eq. (fT2|) . The connection with the 
renormalization group approach arises from the coarse- 
graining assumption after Eq. (fT0|) . where we asserted 
that the rapid oscillation averages to zero. 

Nonlinear Elasticity:- The goal of this section is to derive 
the free energy as a function of the strain tensor, when 
the PFC crystal is deformed under a general deformation 



(16) and, 



where F mn is the deformation gradient. Einstein's sum- 
mation convention is used throughout, except for the in- 
dex j in kj, Aj and Lj. In general, the deformation 
gradient can be written as[l3j 



Fmn — RmpUpn 



(17) 



where R mp is a pure rotation matrix and U pn is a positive- 
definite, pure deformation matrix. Since our system is 
rotationally covariant, we expect that the free energy 
should only depend on the function U T U, where U T is 
the transpose of the matrix U. 

Under the deformation, Eq. (|16p . the complex ampli- 
tudes transform as 



Aj -> Aj = Ae 



(18) 



where we defined D mn = RmkUkn 
that | A,- 1 = A for all j, where A is a constant. kj m 

is the m-th component of the vector kj. Because the 
deformation gradient only enter the complex amplitude 
through its phase, describing local density deformations, 
the only relevant terms in the free energy are the gradient 
terms given by 



E^A'jL*Aj. 

3=1 



(19) 



Other terms in the free energy, Eq. (jlip . only contribute 
when we minimize the free energy with respect to A at 
the end of the calculation. By using Eq. (fT8|) and differ- 
entiating, we obtain 



LjAj ( kj m kjnRmpRnaUpqU a q ~\~ l)^4.j. 



(20) 



Apply Lj again and substitute the result into Eq. JT 
to obtain 



E = A 2 {E 1 - 2E 2 + 3), 



(21) 



where 



E\ I ^ kj rn kj n kj u kj v I F rn q F n q F uw F vw ^ (22) 



and, 



F2 I ^ \ kjm kjn J F m qF n q 



(23) 



The rest of the derivation concerns the evaluation of E\ 
and E 2 . We first evaluate E 2 . By using the definition of 
kj from Eq. (|4j), we obtain 

k\rn, k\n = $my$ny = $ mutiny j (24) 

3 1 V^3 

^2m^2n — ~^^rnx^nx H~ ^rmj ^ (y^mx^ny ~r~ ^nx "h $my)i 

(25) 



3 1 v^3 

(26) 

Combining these three equations we have 
3 3 3 

^ ^ kjtnfejn — ~^3mn(o~mx &my) — ~^mn- (^^) 

Thus, is given by 

£2 = ^p^p, = \u pq U pq = ^Tr[C/ T [/], (28) 

where we used the property of the rotation matrix, that 
RimRjm — RmiRmj = ^[A] and A T are the trace 
and transpose of the matrix A respectively. The evalua- 
tion of Ei is more involved. We note that 

k\mk\nk\uk\v — &mnO'nu$uvO~vnyi (29) 

and observe that the term 

k-imkinki-uhiv + k^ m k^ n k^ u k^ v (30) 
is equal to the term 

2x [k2mk2nk2uk2v+ terms with positive coefficients] 

(31) 

By exploring this relation and using the definition of kj , 
Eq. (|4]), we obtain 

3 3 

^ ^ kjmkjnkj u kj v ~^(^xxyy ~t~ ^yyxx ~t~ ^xyxy) 

3 = 1 

3 

~t~ ~^(^xyyx ~t~ ^yxxy H~ ^yxyx) 

9 

+ g {^xxxx + ^-yyyy): (32) 

where we defined £± a bcd = &ma&nb&uJ>vd for convenience. 
By this, we have 

' 3 \ 
K 3 = l J 

~,[RxpRxaRysRyt H~~ FtypRya^xsRxt\ 

3 

"^[RxpRyaRxsRyt H - RypRxaRys^xt\ 

3 

~g \RxpRyaRysRxt "t - RypRxa-Rxs-Ryt] 
9 

~ T:\Rxp RxaRxsRxt H - RypRyaRysRyt\ (33) 
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To evaluate this expression, we note that we can combine 
terms judiciously. For example, by using the property of 
the rotation matrix, we obtain 

RxpRxaRysRyt "h RxpRxaRxsRxt — RxpRxa&st i (*^) 

and, 

RypRyaRxsRxt 4" RypRyaRysRyt — RypRya^ st • (^^) 

The sum of Eq. (JMJ) and ([35]) then give 

RxpRxaRysRyt RxpRxaRxsRxt 

H~ RypRyaRxsRxt ~t~ RypRyaRysRyt — ^st^ap (^6) 

Repeating this for all the terms in Eq. |[55J) , we obtain 

\ 

^ ^ kjmkjnkju^jv I RmpR-naRusRvt 



[<5ap<5st + S sp S a t + S as Spt] 



(37) 



Substituting into Eq. ([22]). and then into Eq. (f21]) . we 
obtain = 3^4 2 A, where 

A = i{[Tr([/ T J7)] 2 + 2lr(U T UU T U)} - Jr(U T U) + 1. 

(38) 

By using the relation U^ lp U pn — 8 mn + 2u mn , where 
Uij = (diUj + djUi + (dkUi)(dkUj))/2 is the nonlinear 
strain tensor, we obtain 

A = [\ u lx + \ u ly + u ly + u lx + UxxUyyj . (39) 



Finally, we substitute back into the free energy, Eq. 
and minimize the whole expression with respect to 
A to obtain 



A(A) 



-P0±o\/ 9 P0 



i5(r-A) , 



(40) 



which gives the free energy density as 



45 



/(A) = fA 4 (A) + l2p Q A 3 (A) - 3(r - A)A 2 (A). (41) 



This formula completely defines the elastic properties of 
the PFC model, and provides a starting point for con- 
ventional continuum mechanical applications of nonlin- 
ear elasticity theory. 

Nonlinear elastic theories are coarse-grained models 
up to the level of the continuum, and so do not explic- 
itly include defect structures; our results are most use- 
ful for understanding large deformation behaviour, twin- 
ning and phase transitions in ordered materials. How- 
ever, for plastic deformations, dislocations need to be 
included, usually by postulating a free energy with con- 
tributions from nonlinear elasticity (such as Eq. J41[) ). 
vacancies and dislocations [13, While this can yield 
useful insights, the most suitable level of description for 
probing the multi-scale phenomena accompanying plas- 
tic flow remains the PFC equations taking into account 
vacancies fl9| or their rotationally-covariant renormalized 
counterparts (if only dislocations are present). We will 
report on this approach in a future publication [20l| . 
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